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While all materials reduce their intrinsic volume under hydrostatic (uniform) compression, a select few actually expand along 
one or more directions during this process of densification. As rare as it is counterintuitive, such “negative compressibility” 
behaviour has application in the design of pressure sensors, artificial muscles and actuators. The recent discovery of surprisingly 
strong and persistent negative compressibility effects in a variety of new families of materials has ignited the field. Here we 
review the phenomenology of negative compressibility in this context of materials diversity, placing particular emphasis on the 
common structural motifs that recur amongst known examples. Our goal is to present a mechanistic understanding of negative 
compressibility that will help inform a clear strategy for future materials design. 


1 Introduction 

Negative linear compressibility (NLC) is the bizarre materials 
property whereby a system expands along one direction when 
compressed uniformly. 1 2 Bizarre, because our intuition is that 
materials should shrink when squeezed—an intuition that is 
(rightly) grounded in the thermodynamic requirement that vol¬ 
ume be reduced at increased pressure. 3 4 Yet NLC does not vi¬ 
olate thermodynamics: it simply arises whenever volume re¬ 
duction can be coupled to linear expansion [Fig. [TJ. In the 
benchmark review of NLC—now 17 years old—Baughman 
explains how the phenomenon might eventually be applied 
in a variety of ways, including the development of artificial 
muscles and amplification of piezoelectric response for next- 
generation sensors and actuators. 2 Until recently, there has 
been relatively little hope of identifying suitable candidates 
for these applications. The most significant challenges have 
been the apparent rarity of NLC (Ref. |2 reports it to occur 
in only 13 known materials) and the extreme weakness of the 
NLC effects exhibited by these materials. 

Over the past few years, the field has changed in two 
important respects. The first is that materials have now 
been discovered that exhibit orders-of-magnitude stronger 
NLC effects than the “classical” NLC materials reviewed by 
Baughman. 2 The second advance—which likely reflects the 
improved accessibility of variable-pressure crystallographic 
measurements—is that NLC has now been found to occur in a 
much greater diversity of materials, ranging from dense inor¬ 
ganic oxides^ and fluorides® to metal-organic frameworks®^ 
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Fig. 1 Mechanical responses to hydrostatic pressure: (a) positive 
compressibility—contraction in all directions; (b) negative linear 
compressibility—linear expansion in one direction; (c) inflation 
associated with incorporation of the pres sure-transmitting media 
(blue circles) within the material interior. The system volume 
(represented here by the solid red area) is reduced in all cases. 


and even molecular solids. 9 Consequently we felt it timely 
to review the phenomenon of NLC once again, placing par¬ 
ticular emphasis on the common underlying geometric motifs 
responsible for NLC in the various materials—both old and 
new—and in doing so to help inform future materials design. 

Our review is organised as follows. We begin with an 
overview of the theoretical and experimental approaches to 
understanding, measuring and comparing NLC responses. As 
part of this discussion we introduce the new measure of “com¬ 
pressibility capacity”, Xk- This will play a role in allowing us 
to compare the NLC responses of very different materials. The 
bulk of the review concerns the NLC behaviour of known ma¬ 
terials, grouped according to the microscopic mechanism re¬ 
sponsible for NLC. The summary with which our review con¬ 
cludes aims to collate succinctly the various data presented, 
making particular use of the Ashby plot approach. We also 
discuss the design criteria for different applications of NLC 
materials and summarise the various directions in which we 
expect the field to develop over the coming years. 
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2 Compressibility: Theory and Measurement 


In the simplest terms, the compressibility of a material de¬ 
scribes the relative rate of collapse of its linear dimensions 
with respect to pressure, measured or calculated at constant 
temperature : 3 4 


Ki = - 



( 1 ) 


The minus sign means that positive compressibilities corre¬ 
spond to length reduction under increasing pressure. Con¬ 
ventional engineering materials such as steel and concrete 
contract by ^0.5% in every direction for each GPa of ap¬ 
plied pressure, corresponding to a linear compressibility K ~ 
5TPa -1 . 10 Compressibility magnitudes usually reflect bond 
strengths, and so softer materials such as polymers and foams 
exhibit much larger values; for example, the linear compress¬ 
ibility of polystyrene is K ~ lOOTPa -1 . 11 

Crystalline materials will in general have different com¬ 
pressibilities in different directions. For example, a layered 
material will usually be more compressible along the stacking 
axis than it is along a perpendicular direction . 12 This direc¬ 
tional dependence can be relatively complex, especially when 
the crystal symmetry is low. We proceed to introduce the the¬ 
ory of compressibility in its most general form before explain¬ 
ing how the situation can be simplified as symmetry increases. 
Our starting point is the formal definition of compressibility as 
a rank -2 tensor : 3 


K = — 
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( 2 ) 


Here the £^- are functions of hydrostatic pressure p and repre¬ 
sent the pressure-induced strain experienced by axis j along 
axis i. The eigenvectors of Eq. describe an orthogonal 
coordinate system that brings K into diagonal form. These 
vectors are the so-called “principal axes” of compressibility 
(sometimes labelled x \, X 2 , X 3 ) which can be interpreted as the 
crystal directions along which hydrostatic compression does 
not lead to any shear component. The eigenvalues of K, which 
we term ^ 4 ,£ 2 ,£ 3 , correspond to the compressibilities along 
these principal axes and are the unique descriptors of linear 
compressibility for any crystalline material. The formal re¬ 
quirement for NLC is that at least one of the A) is negative. 

Defined in this way, the linear principal compressibilities 
are directly related to the volume compressibility, and in turn 
to the bulk modulus: 



B = Ky l 


K l +K 2 + K 3 , (3) 


K l +K 2 + K 3 ’ 


Because the volume compressibility must be positive, any sys¬ 
tem for which one of the linear compressibilities exceeds the 
bulk compressibility (i.e., K[ > Ky = B~ l ) must exhibit NLC. 
This is the type of approach to identifying NLC materials em¬ 
ployed in Ref. 0 

(A brief aside—Conventions vary in terms of the symbols 
used to denote these various elastic parameters. Compressibil¬ 
ities are denoted by some using the symbol (3 which is used 
by others to mean the volumetric coefficient of thermal expan¬ 
sion, 13 and by perhaps very many more to mean one of the unit 
cell angles. Likewise the bulk modulus is denoted by K within 
much of the mineralogical literature, despite this symbol as¬ 
suming the inverse meaning of compressibility when used in a 
physics text.® In this review we adopt the conventions of the 
condensed matter physics community— i.e. K for compress¬ 
ibility and B for bulk modulus—which we feel are the least 
likely to cause confusion.) 

2.1 Compressibilities from variable-pressure crystallo¬ 
graphic measurements 

As mentioned above, the tensor algebra associated with com¬ 
pressibility determination is simplified enormously by consid¬ 
eration of crystal symmetry. For systems of orthorhombic 
crystal symmetry or higher, the principal axes coincide with 
the crystal axes. This means that the lattice parameter com¬ 
pressibilities 
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( 6 ) 

(7) 


which can be determined using variable-pressure crystallo¬ 
graphic measurements, give directly the principal axis com¬ 
pressibilities. In other words, the Ki reflect the relative rate of 
change of the lattice parameters with respect to pressure, and 
NLC materials can be identified as those for which at least one 
lattice parameter increases under hydrostatic pressure. 

Unfortunately this equivalence between lattice and princi¬ 
pal axis compressibilities does not hold for systems with mon¬ 
oclinic or triclinic crystal symmetries; lattice parameter com¬ 
pressibilities can have very little direct physical meaning in 
these cases . 15 In particular a negative value of one or more lat¬ 
tice parameter compressibilities would no longer imply NLC 
because the principal axis compressibilities may nonetheless 
remain positive. For such situations, there are software pack¬ 
ages that facilitate the conversion from lattice parameter to 
principal axis compressibilities: PASCal and EoSfit are two 
examples . 15 16 
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One further complication in converting variable-pressure 
lattice parameters to linear compressibilities is the tendency 
for lattice parameters to depend non-linearly on pressure; i.e. 
the K[ are themselves pressure-dependent. (Indeed if any 
K[ were truly constant then there would exist a finite pres¬ 
sure at which the corresponding material length would van¬ 
ish: p cr it = 1 /K[). There is no thermodynamic requirement for 
the K[ to depend on pressure in any particular way; this sit¬ 
uation contrasts that of the pressure-dependence of the crys¬ 
tal volume, which is often int erpr eted in terms of the Birch- 
Murnaghan equations of state. 17 18 Instead the lattice param¬ 
eters are usually fitted to some empirical parameterisation of 
choice. In the simplest case this would be the linear relation¬ 
ship 

£(p)=to[l-K £ p\, (8) 

where £o represents the length at zero pressure. The nonlinear¬ 
ity left unaccounted for by this simple parameterisation can be 
included via higher-order polynomial expansions: 


(■(p) = 4 


1 + £ Clip' 
r=i 


(9) 


Some authors then identify the value — oq with the linear com¬ 
pressibility Kf \ however, this tends to overestimate the com¬ 
pressibility in cases where there is a strong pressure depen¬ 
dence. 19 2 din our own work, we have found that the alterna¬ 
tive parameterisation 


t(p) = 4 + A(p-p c ) v (10) 


actually captures better the pressure-dependence of lattice pa¬ 
rameters for most systems. The corresponding compressibili¬ 
ties are determined straightforwardly from the pressure deriva¬ 
tive of Eq. 


K (p) = ~j^kv( p -p c ) v (11) 

The value of K(p) determined in this way diverges at p = p c , 
since V < 1; often this divergence has physical significance 
in terms of the elastic instability at a pressure-induced phase 
transition. Figure [2] compares the linear compressibilities 
determined in these different ways for some representative 
variable-pressure lattice parameter data. 

Whichever parameterisation is used, accurate determination 
of the Ki generally relies on access to a relatively large number 
of lattice parameter measurements over the pressure interval 
of interest. Experimental claims of NLC are sometimes made 
on the basis of just two measurements; 23 however a general 
rule of thumb is that 10 measurements are needed for accu¬ 
rate compressibility determination with little improvement for 
more than 20 measurements. 24 



pressure (GPa) 

Fig. 2 Extraction of linear compressibilities from variable-pressure 
lattice parameter measurements. The top panel shows representative 
lattice parameter data for a smoothly-compressible material—in this 
case measured for KMn[Ag(CN) 2 l 3 - 22 The three common 
approaches to fitting lattice parameter data represented by 
Eqs. ([8l-fT0t give the fits to data shown in red, blue, and green, 
respectively. The corresponding compressibilities are shown in the 
lower panel. All three methods obtain comparable average 
compressibility values over the entire pressure range for which data 
are fitted. The reduction in magnitude of K at higher pressures is 
unaccounted for by the linear fitting method, and is treated slightly 
differently by the two non-linear fits. The extrapolated values of K 
at p = 0 are very different in all three cases. 


2.2 Compressibilities from elastic compliances 

Because lattice compressibilities are one aspect of the more 
general elastic behaviour of materials, determination of the 
elastic stiffness tensor C provides an alternative means of 
characterising NLC. 4 The relationship between the Kf and C 
is most straightforwardly established by considering the ele¬ 
ments of the elastic compliance tensor S (the inverse of C), 
which relate strains to applied stresses Gif 

3 

£/j = SijklGkl- ( 12 ) 

k,l= 1 
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Here all terms are defined in the limit of infinitesimal strain. 
Eq.© is essentially a form of Hooke’s law that relates dis¬ 
placement (oc strain) to mechanical force (oc stress). Hydro¬ 
static compression is the specific situation where axial stresses 
are all equal to the applied pressure and shear stresses are ab¬ 
sent; that is, 

o u = p Vz, (13) 

Qy = 0 (14) 

By design, the principal axes are those for which axial com¬ 
pression does not induce any shear strain, and so Eq. © re¬ 
duces to 

3 

£ii = ~P ^ Sik ( 15 ) 

k= 1 

when expressed in the principal axis coordinate system (not¬ 
ing that we have now switched to Voigt notation so that S# = 
Snkk)- Substitution into Eq. ([2]) gives the simple relationship 

3 

K, = £ Sij. (16) 

7=1 


The tensor S is obtained by inversion and the K L calculated ac¬ 
cording to Eq. ( p~6l ). Such an approach is implemented in, for 
example, the CRYSTAL09 code,® 34 and has been applied to 
the exploration of negative compressibility of simple inorgan¬ 
ics] 35 36 zeolites, 32 1 and metal-organic frameworks alike. 37 38 
Once again, the compressibilities obtained in this way rep¬ 
resent the zero-pressure limit and as such can be vastly more 
extreme than those determined across finite pressure ranges J 37 ^ 
This link between linear compressibilities and the more 
general elastic properties of materials means that the obser¬ 
vation of NLC is often diagnostic of other anomalous me¬ 
chanical responses. Examples include extreme mechanical 
anisotropy/ 37 unusual valu es of the Poisson’s ratio, 19 negative 
thermal expansion (NTE)J 12 39 and a propensity towards dy¬ 
namic instabilities. 19 -l Importantly, this correspondence works 
both ways. Severe structural anisotropy (such as arising from 
the preferred orientations of molecules in a particular packing 
arrangement or the symmetry of a given framework topology), 
or the observation of uniaxial or biaxial NTE are increasingly 
frequently found to be strong predictors of NLC behaviour. 8 40 
We will come to explore this correspondence in more detail 
below. 


that often appears in texts on the subject. 4 Consequently, de¬ 
termination of the elastic stiffness tensor C (either experimen¬ 
tally or computationally) can also yield the linear compress¬ 
ibilities Ki via the compliance tensor S = C -1 . 

Experimental techniques capable of probing the tensor C 
include resonant ultrasound spectroscopy (RUS)J 25 Brillouin 
scattering, 2627 inelastic neutron scattering, 28 nanoindenta¬ 
tion 29 and shear-wave velocity 30 measurements. In most 
cases these measurements are indirectly sensitive to a subset 
of the elements, or combinations of elements, of C. Conse¬ 
quently interpretation of the experimental data is often carried 
out via a parameterised lattice-dynamical model (containing 
fewer free parameters than C itself), from which C is subse¬ 
quently calculable using software packages such as GULP.cn] 
The single most important distinction between experimen¬ 
tal compressibilities determined via elastic compliances and 
those obtained from variable-pressure crystallographic mea¬ 
surements is that the former correspond to values obtained in 
the limit of zero applied pressure and hence are usually much 
larger in magnitude. 

Ab initio methods also allow determination of uniaxial com¬ 
pressibilities in an analogous way. Starting from the optimised 
geometry for a given crystal, sampling of every possible pair¬ 
wise combination of strains £ij,£jd allows the Qjki to be evalu¬ 
ated directly from the corresponding second derivatives of the 
lattice energy: 32 


Qjki = - 


V 


d 2 E 

d £i j d £ki 


(17) 


2.3 Compressibility capacity 

Because of the pressure dependence of the K t it is important 
that compressibility values are quoted with reference to the 
pressure range over which they have been determined. This 
pressure range may correspond to the entire stability field of 
the phase in question, may be imposed by experimental con¬ 
straints, or may be of relevance to a specific industrial pro¬ 
cess (by way of example, most machining processes subject 
materials to pressures of ca 1-2 GPa®). One obvious limita¬ 
tion in comparing the NLC behaviour of different materials in 
terms of the magnitudes of K t alone is that these values may 
be determined over very different pressure ranges for different 
materials. 

A metric that we have found useful in comparing the degree 
of NLC behaviour for different materials is what we term the 
compressibility capacity: 

Xk = ~ f PC K(p) dp, (18) 

Jo 

where p c denotes the maximum pressure for which NLC is 
observed. The value of Xk is a dimensionless quantity that 
takes into account both the magnitude of NLC and the pres¬ 
sure range over which NLC occurs. It simply represents the 
maximum total fraction by which a crystal can expand under 
application of hydrostatic pressure. We will report Xk values 
for the various NLC materials reviewed in this article as and 
when each is introduced; anticipating these results, we find 
that Xk 1% for the majority of known NLC systems, with 


4 | Journal Name, 2010, [volfl - [i~8| 


This journal is ©The Royal Society of Chemistry [year] 









values closer to 10% only observed in the most exceptional 
cases. 


3 NLC Materials 

What follows is an overview of the various materials known to 
exhibit NLC. We have grouped these into four classes accord¬ 
ing to the microscopic mechanism likely to be responsible for 
NLC in each case: (i) those compounds for which NLC arises 
as a consequence of proper, improper, or quasi-ferroelastic 
phase transitions; (ii) network solids for which NLC is driven 
by correlated polyhedral tilts; (iii) helical systems; and (iv) 
framework materials with wine-rack, honeycomb, or related 
topologies, where NLC arises from framework hinging. In this 
way, we hope to summarise not only the properties of known 
NLC materials but also the mechanisms that can give rise to 
the phenomenon itself. 


3.1 Ferroelastics 


Ferroelasticity describes the emergence of spontaneous strain 
in a symmetry-breaking phase transition. 42 The phenomenon 
can be considered the mechanical equivalent of ferroelectricity 
or ferromagnetism, where it is the spontaneous strain—rather 
than polarisation or magnetisation—that behaves as the ferroic 
order parameter. A simple example is the square —> rhombic 
transition illustrated in Fig. [3j for which the ferroelastic or¬ 
der parameter is a measure of the distortion away from square 
symmetry:®® 


0 = 


b — a 
a + b 


(19) 


The value of (j) is zero in the ‘paraelastic’ high-symmetry 
phase and non-zero in the ferroelastic low-symmetry phase 
(square and rhombic geometries in Fig. [3] respectively). One 
result of Landau theory is that any transition for which strain is 
the primary order parameter must be second-order in nature. 43 
So, there being no volume discontinuity across these so-called 
‘proper ferroelastic transitions’, the emergence of spontaneous 
strain (i.e. (j) ^ 0 or, equivalently, a ^ b in Eq. (19] )) requires 




Fig. 3 The paraelastic (left) and ferroelastic (right) states of the 
square lattice. The ferroelastic order parameter 0 = (b — a)/(a + b) 
is a measure of the extent of symmetry-breaking distortion. 




HCP (NiAs-type) 



Fig. 4 Ferroelastic NLC mechanism in rutiles, (a) The ferroelastic 
instability of the rutile structure type corresponds to a progressive 
distortion towards a defect NiAs structure and involves expansion in 
the vertical direction. 46 ^ (b) Pressure-dependence of the lattice 
parameters of TeC> 2 , showing the ferroelastic transition at 
- 0.81 GPaJElThe high -pressure CaCL-structured phase exhibits 
NLC along b. (c) Evolution of the octahedral tilt parameter ft>. 4 ^(d) 
Softening of the effective elastic constant C' at the ferroelastic 
transition. 30 


at least one axis to increase in length on symmetry lower¬ 
ing. This means that we can expect NLC in any system that 
supports a pressure-induced proper ferroelastic transition: the 
lower-symmetry phase stabilised at pressures above the tran¬ 
sition must expand along at least one crystallographic axis for 
at least some finite pressure interval. 


Rutiles 

A number of simple binary inorganic solids with the TiC >2 
rutile structure exhibit NLC via precisely this mechanism.® 
The well-known tetragonal crystal structure of rutile con¬ 
sists of columns of edge-sharing octahedra which are in turn 
connected at their comers [Fig. [4ja )]. This structure has a 
ferroelastic instability associated with correlated rotations of 
neighbouring columns of octahedra. 30 48 Activation of this 
tilt system lowers the crystal symmetry from tetragonal to 
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Table 1 Compressibilities of materials for which NLC arises as a consequence of ferroelastic or related phase transitions. 



Method 

K\ (TPa" 1 ) 

k 2 (TPa- 1 ) 

K 3 (TPa" 1 ) 

B 0 (GPa)“ 

Range (GPa) 

Xk ( % ) 

Ref(s). 

Te0 2 

a 

-5.1(6) 

2.1(7) 

18.4(6) 

52(4) 

0.9-3.25 

1.20(14) 

MMWWWWM 

NiF 2 

a 

-0.48 

0.61 

4.41 

222 

1.8-3.2 

0.067 

mn 

/3-Mn0 2 

a 

-0.16(7) 

0.269(17) 

1.82(10) 

328(18) 

0.3-29.3 

0.46(20) 

M 

MgF 2 

a 

-1.3(3) 

2.51(2) 

8.05(10) 

68(13) 

9.1-10.4 

0.17(4) 

52 53 

Pb0 2 -I' 

a 

-1.82(15) 

1.76(7) 

3.40(3) 

167(18) 

3.8-6.1 

0.42(3) 

ED 

Ge0 2 

a 

-0.137(12) 

0.827(14) 

2.02(13) 

— 

28-36 

0.110(10) 

\55\ 

Sn0 2 

a 

-0.185 

0.394 

2.345 

204(6) 

11.8-21 

0.170 

145156 

Zn(CN) 2 TI 

b 

-2.08 

11.2 

16.8 

— 

1.52-5 

0.724 

|57|58 

Pb 3 (P0 4 ) 2 

a 

-4.3(4) 

5.0(2) 

20.98(18) 

38.7(5) 

0-1.59 

0.68(6) 

E9HEQ 

InS 

a 

-2.41(13) 

2.8(4) 

15.3(22) 

33.2(18) 

0-4.3 

1.04(6) 

mm 

Sillimanite 

a,c 

-3.30 

1.45 

10.8 

112.50 

29.9-37.5 

2.51 

m 

PtS 

b,c 

-0.47 - -0.92 

3.3 

3.3 

167.5 

0-10 

0.47-0.92 

36)66 


a Calculated from PASCal 15 from reported variable-pressure lattice parameters. b As reported. c From DFT calculation. 


orthorhombic, resulting in the (equally well-known) CaCl 2 
structure [Fig.|4ja)]. In order to conserve volume, the lattice 
is forced to expand along one of the two directions perpendic¬ 
ular to the column axisM 7 67 This transition can be viewed as 
a progression towards the defect-NiAs-type structure of FeS 2 
marcasite [Fig. |4ja)] and is thought to be driven largely by 
considerations of anion packing efficiency.® 

A representative example of a pressure-induced ferroelas¬ 
tic transition—and hence NLC—arising from this ty pe of in¬ 
stability is given by the mineral paratellurite (Te 02). 5 47 49 50 
Variable-pressure lattice parameter measurements for this ma¬ 
terial reveal a ferroelastic transition at a critical pressure p c ~ 
0.9 GPa, followed by NLC along the b crystal axis of the re¬ 
sulting high-pressure phase [Fig.|4jb)]. The proposed mech¬ 
anism of correlated octahedral tilts can be validated by deter¬ 
mining the pressure-dependence of the octahedral tilt angle 
CO , which also behaves as an order parameter for the transi¬ 
tion [Fig.|4jc)].& 8 That the transition is truly strain-driven can 
be deduced from the elastic behaviour near p c : the effective 
elastic constant C — \{C\\— Cn) governs the relevant shear 
mode velo city, a nd can be seen to vanish at the transition point 
[Fig.^d)].™ 

TeC >2 is not an isolated example of this behaviour: a num¬ 
ber of rutile-structured dioxides and difluorides exhibit NLC 
via essentially the same mechanism.® The magnitude of 
NLC for these different materials can be determined from the 
various crystallographic measurements reported in the litera¬ 
ture. In the case of TeC> 2 , the data of Ref. l50i give £nlc = 
—5.1(6) TPa -1 over the pressure range 0.9-3.25 GPa, corre¬ 
sponding to a compressibility capacity of %k = 1.20(14)%. 
Table [I] compares these values for the various rutiles for 
which a high-pressure ferroelastic phase transition has been 
observed. What emerges is that there is a general corre¬ 
spondence between cation radius and extent of NLC, such 
that TeC> 2 —which contains the largest® of the cations—also 



Fig. 5 Relationship between magnitude of NLC and cation radius 
for rutile-structured dioxides and difluorides. Radii taken from Ref. 
l69l The line is a guide to the eye. 


shows the most extreme NLC response [Fig.[5). We will show 
by comparison with other families that a compressibility ca¬ 
pacity Xk — 1% is not itself particularly extreme . So while it 
is the case that a ferroelastic instability mechanism for NLC 
may give rise to relatively general behaviour, it would seem 
that—at least in the case of rutiles—the mechanism is unlikely 
to give rise to especially large NLC responses. One additional 
complication is that the NLC effect is of course only observed 
in the non-ambient phase, so future research within this fam¬ 
ily might likely concentrate on lowering p c (perhaps even to 
negative pressures; i.e., studying CaCl 2 -structured materials) 
and/or varying cation/anion radii so as to allow the largest 
possible spontaneous strains to emerge within the ferroelastic 
phase. 

Zinc cyanide: an improper ferroelastic 

Even materials for which hydrostatic pressure induces a so- 
called ‘improper’ ferroelastic phase transition can exhibit 
NLC within the high-pressure phase. The label ‘impoper’ 
simply means that the spontaneous strain is no longer the pri- 
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mary order parameter responsible for driving the phase transi¬ 
tion. 42 Instead ferroelastic strain develops via coupling to an 
alternate, dominant, symmetry-breaking mechanism, such as 
a phonon instability . 70 The Landau conditions change in these 
situations such that volume discontinuities may be observed, 
and so there is no strict guarantee of NLC. Nevertheless if the 
volume collapse on symmetry lowering is small with respect 
to the subsequent evolution of ferroelastic spontaneous strain, 
then NLC may indeed arise. 


This is precisely the mechanism that appears to be respon¬ 
sible for NLC in the high-pressure phase of t he molecular 
framework material zinc cyanide, Zn(CN) 2. 57 58 71 Under am¬ 
bient conditions Zn(CN )2 adopts a cubic structure in which 
Zn 2+ cations are tetrahedrally coordinated by four cyanide an¬ 
ions, each of which in turn connects two Zn centres to give 
a three-dinemensional framework with the diamond-like anti¬ 
cuprite topology . 72 Both computational studies and inelastic 
neutron scattering measurements point to the existence a large 
family of soft phonon modes, many of which are also im¬ 
plicated in the strong NTE behaviour observed experimen¬ 
tally . 71 73 74 Hydrostatic compression to 1.52 GPa results in a 
first-order displacive phase transition to a denser orthorhom¬ 
bic structure, which is consistent with the condensation of 
at least one zone-boundary soft mode . 57 58 The correspond¬ 
ing atomic displacements involve correlated rotations of con¬ 
nected pairs of Zn(C/N )4 tetrahedra, resulting in coupled ex¬ 
pansion/contraction of the crystal lattice perpendicular to the 
rotation axis in much the same way as described for the rutile- 
structured systems discussed above [Fig. [ 6 ^a )] . 57 NLC per¬ 
sists in the high-pressure Zn(CN) 2 -II phase from the I/II tran¬ 
sition at 1.52 GPa up to ^5 GPa and is reasonably strong 
over this entire pressure range: £nlc = — 11(3) TPa -1 and 
Xk = 3.8(10)%. 


In competition with any ferroelastic NLC mechanism is the 
tendency for framework buckling, which favours PLC and 
becomes increasingly important at higher pressures . 39 In the 
specific case of Zn(CN) 2 -II this buckling involves a systematic 
distortion of Zn(C/N )4 polyhedra and bending of Zn-C-N-Zn 
linkages to allow additional neighbouring cyanide ions within 
the originally-tetrahedral Zn coordination sphere . 57 The Zn 
coordination number progressively increases from four to six, 
and the resulting volume reduction becomes increasingly sig¬ 
nificant with respect to the NLC effect of correlated tilts, 
such that the crystal axis along which NLC is initially ob¬ 
served begins to contract for pressures higher than 5 GPa. 
Such crossover between NLC and PLC behaviour might be 
expected to be a reasonably general phenomenon, and has cer¬ 
tainly be noticed previously in e.g. a variable-pressure study of 
the metal-organic framework silver(I) methylimidazolate . 39 



Fig. 6 NLC from improper and reverse ferroelastic transitions, (a) 
The cubic —» orthorhombic improper ferroelastic transition in 
Zn(CN )2 involves correlated rotations of column pairs of Zn(C/N )4 
tetrahedra (shaded in green), resulting in NLC along the a axis of 
the daughter cell . 57 (b) The ferroelastic state of Pb 3 (P 04)2 is denser 
than the paraelastic parent, and so application of pressure induces a 
reverse ferroelastic transition. The c parameter of the ambient phase 
(the length of which is normalised here d = c/yj 3 for comparison) 
expands as the transition is approached.® Pb atoms are shown as 
large black spheres and PO 4 units as filled tetrahedra. 


Reverse ferroelastics: lead phosphate and indium sulfide 

Lead phosphate, Pb 3 (PC> 4 ) 2 , is the unusual example of a 
material that exhibits a ‘reverse’ ferroelastic phase transi¬ 
tion on compression: the low-symmetry monoclinic struc¬ 
ture converts to a higher-symmetry, denser, rhombohedral 
phase at a hydrostatic pressure of approximately 1.6 GPa . 75 76 
Because ferroelastic transitions couple expansion and con¬ 
traction along orthogonal axes, NLC is actually expected ir¬ 
respective of the direction in which the transition between 
high symmetry and low symmetry states is traversed. In¬ 
deed Pb 3 (P 04)2 shows NLC throughout the entire stability 
field of the ferroelastic phase with K^lc = —4.3(4) TPa -1 and 
Xk = 0.68( 6) %.® 77 The mechanism responsible for NLC— 
which resembles that observed on heating Pb 3 (PC> 4)2 at am¬ 
bient pressure—again involves correlated polyhedral tilts (in 
this of the PO 4 tetrahedra) which couple to off-centering of 
the Pb 2+ cations [Fig. [bjb)]. Doping with Ba 2+ predictably 
favours the rhombohedral (cation-centred) state, lowering at 
once both transition pressure and temperature while preserv- 
ing NLC.EEMmi] 

A less-well characterised example of a possible reverse 
ferroelastic transition is that of indium(II) sulfide, InS [= 
(In 2 ) 4 + (S 2 _ ) 2 ].^® Certainly its ambient phase (orthorhom- 
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bic Pmnn symmetry) shows NLC over the pressure range 
0-4.3 GPa (£nlc = — 2.41(13) TPa -1 ), and this behaviour is 
qualitatively understandable in terms of the straightening of 
S-In-In-S “dumbells”. The confusion lies in the relationship 
of this response to the high-pressure InS-II phase that forms at 
7.5 GPa. This high-pressure phase was originally reported to 
have the tetragonal Hg 2 Cl 2 structure (i.e. with linear S-In-In- 
S units), which is of exactly the right symmetry to be consid¬ 
ered the paraelectric parent of the ambient phase. 62 However, 
a more recent in situ single crystal structure determination re¬ 
ports InS-II to adopt a distorted monoclinic structure in which 
the same S-In-In-S units remain buckled.® So while there 
is no ambiguity regarding the NLC effect itself, and there is 
an implied relationship to the existence and nature of a high- 
pressure phase transition, a detailed mechanistic understand¬ 
ing of this relationship would demand further experimental 
characterisation of the pressure-dependent behaviour of this 
system. 


NLC from ferroelastic-like phase transitions 

Other types of phase transition—which are not strictly fer- 
roelastic, but which share mechanistic similarities—may also 
give rise to NLC. We illustrate this point with two final ex¬ 
amples. The first involves a geologically-relevant isosym- 
metric transition in the sillimanite polymorph of Al 2 Si 05. 65 
The ambient phase of this framework structure contains four- 
coordinate Si atoms. 78 Lattice dynamical calculations using 
well-established interaction potentials optimised for alumi¬ 
nosilicates 79 indicate that the framework undergoes a corre¬ 
lated reorganisation between 30 and 40 GPa that increases the 
Si coordination number from four to five. 65 All original bond¬ 
ing connectivity is preserved in this process (only new bonds 
are formed), and there is no change in crystal symmetry; such 
a transition is necessarily first-order in natureJ 80 81 As a re¬ 
sult of the increased Si coordination number, the high-pressure 
phase is denser than the low-pressure phase, but the new con¬ 
nectivity results in an increase in length along the c axis of its 
orthorhombic Pnma cell. As the transition is approached with 
increasing hydrostatic pressure, the gradual conversion of one 
phase to the other results in an NLC effect [Fig.[7Ja)]. 

Our second example is the transformation of PtS from its 
ambient-pressure structure (tetragonal P42/nmm) to a high- 
pressure phase with the so-called PdS structure (tetragonal 
P42m) at pressures of between 2.5 and 3.0 GPa.® From an 
experimental viewpoint, very little is understood regarding 
this transition. What is known is that the high-pressure phase 
is metastable under ambient conditions and is expanded by 
^8% along the tetragonal axis and compressed ^9% along 
the two perpendicular axes relative to the thermodynamic 
phase.® First principles calculations suggest that the defor¬ 
mation mechanism of the ambient phase under compression 





Fig. 7 (a) The first-order isosymmetric transition in sillimanite 
involves a discontinuous increase in Si coordination number from 
four to five (SiO„ polyhedra shown in blue). The corresponding 
lattice rearrangement involves an expansion along the c axis.^ 5 (b) 
At hydrostatic pressures of ~3-4 GPa, PtS undergoes a 
reconstructive transition to the denser PdS structure type. 36 ® Both 
phases have tetragonal symmetry but the latter expanded along c 
relative to the former. Quantum mechanical calculations show that, 
as the transition is approached on increasing pressure, the ambient 
phase is expected to expand along the same axis.® 


resembles a progression towards this high-pressure structure, 
which means that coupled PLC/NLC is expected, with NLC 
occurring along the tetragonal axis [Fig. [7^b)]. 36 The magni¬ 
tude of NLC obtained in these calculations turns out to de¬ 
pend on the particular functionals used: the authors of Ref. 
[36] find K nlc = -0.47 TPa” 1 for LDA and —0.92 TPa -1 for 
GGA (calculated over the pressure range 0-10 GPa in both 
instances). The large expansion between the two phases ob¬ 
served experimentally (which places an upper limit on the 
true value of %k) suggests that the real value of K^lc may 
be higher than these calculated values, and clearly additional 
experimental characterisation would play a valuable role in 
understanding better the intriguing NLC behaviour of this sys¬ 
tem. 

So the structural changes that occur near phase transi¬ 
tions can give rise to NLC for a variety of different struc¬ 
tural families and different transition mechanisms. NLC is 
all but guaranteed in the case of proper ferroelastic transi¬ 
tions, given that the development of spontaneous strain with¬ 
out a volume discontinuity requires expansion of the crystal 
lattice along at least one direction. Ambient-pressure variable- 
temperature studies may provide a useful method of identify¬ 
ing likely NLC candidates for further study, by virtue of the 
empirical observation that ferroelastic instabilities observed 
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on cooling are also often observed under hydrostatic pres¬ 
sure. In this respect the recent discovery that some ferroelas- 
tic metal-organic frameworks (MOFs) can develop extremely 
large spontaneous strains on cooling suggests that equally ex¬ 
treme NLC may also be discovered in the very same sys¬ 
tems.^ 


3.2 Tilting networks 

The concept that correlated polyhedral tilts might give rise to 
NLC even in the absence of a phase transition is by no means 
new—and is actually more likely to yield practically useful 
systems since NLC is then an intrinsic property of the am¬ 
bient phase. In the nonlinear optic (NLO) material BiB 3 C> 6 , 
for example, the dominant deformation mechanism under hy¬ 
drostatic pressure involves correlated tilts of BO 3 units which 
act to hinge the connected borate framework.ESMIxhis mech¬ 
anism drives a strong contraction in one direction (parallel 
to the highly-compressible Bi 3+ lone pairs) that couples to a 
moderate expansion in a perpendicular direction [Fig.^a )]. 83 
There is no change in crystal symmetry throughout this pro¬ 
cess; in particular the mechanism cannot be thought of as aris¬ 
ing in the vicinity of a ferroelastic phase transition. 

A further point of interest regarding BiB 306 is that it is 
one of the few NLC systems for which compressibilities have 
been determined using two complementary techniques. On the 
one hand, a variable-pressure crystallographic study reported 
^nlc = —6.7(3) TPa -1 over the pressure range 0-5 GPa (i.e., 
Xk = 3 .35(15)%).^ 3 On the other hand, experimental determi¬ 
nation of the elastic tensor based on RUS measurements gave 
A^nlc = —12.5 TPa -1 , which is relevant in the limit p —>• 0 
[Fig. ^b)].® The factor-of-two difference between the two 



Fig. 8 (a) Correlated tilts of BO 3 units polyhedra lead to hinging of 
the borate framework in B^Of, and volume reduction; the network 
expands along a in the process.!- 83 Bi atoms are shown as isolated 
black circles and the BO 2 network in polyhedral representation, (b) 
The compressibility indicatrix^ determined from RUS 
measurements shown in the same orientation as (a): red and blue 
regions indicate, respectively, positive and negative values of the 
linear compressibility.^ 



Fig. 9 (a) Correlated rotations of PO 4 tetrahedra about the b axis 
would cause the CSH 2 PO 4 structure to densify while expanding 
along a direction parallel to c. 87 H atoms are shown as white 
spheres, PO 4 units in polyhedral representation, and Cs atoms 
omitted for clarity.^ (b) The densification mechanism of a generic 
herringbone lattice, which acts as a geometric model for NLC in 
CsH 2 P0 4 . 


measurements of K^lc does not imply an experimental incon¬ 
sistency, but rather reflects the real variation in K with pressure 
as illustrated in Fig. [2] That NLC is relatively strong and per¬ 
sistent in BTB 3 O 6 likely reflects the openness of its framework 
structure and the polarisability of the extra-framework Bi 3+ 
cation. 

Equally strong NLC effects are expected to occur amongst 
other open framework structures based on connected polyhe¬ 
dra. In a recent computational study of 121 siliceous zeolites, 
a total of 16 were identified as NLC candidates on the ba¬ 
sis of their calculated elastic compliance tensors . 32 Amongst 
these frameworks, the strongest NLC behaviour is anticipated 
for the experimentally-realisable GIS zeolite topology , 86 with 
Xnlc = —13.7 TPa -1 in the p —» 0 limit. What is remarkable 
here is the unexpected frequency of NLC: if the phenomenon 
occurs in 13% of a broad family of materials then it is certainly 
less rare than originally envisaged . 32 The pressure range over 
which these zeolites exhibit NLC remains to be determined; 
likewise the effects of A1 substitution, extra-framework coun¬ 
terion inclusion, and solvation on NLC offer additional av¬ 
enues of experimental investigation. 

Correlated tilts of PO 4 tetrahedra have also been implicated 
in the strongly anisotropic elastic response of the widely- 
studied ferroelectric material CSH 2 PO 4 [Fig.|9ja)].S®To the 
best of our knowledge, the only variable-pressure crystallo¬ 
graphic measurements of the ambient phase of this material 
are limited to (i) single-crystal neutron diffraction studies of 
the evolution of diffuse scattering , 91 and (ii) low resolution X- 
ray powder diffraction patterns at 0.29 and 2.89 GPa reported 
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Table 2 Compressibilities of materials for which NLC arises as a consequence of correlated polyhedral tilts. 



T^nlc (TP a A ) 

74>lc (TPa -1 ) 

B (GPa) 

Range (GPa) 

Xk (%) 

Ref. 

BfELOft 

-6.7(3) 

- 


0-5.0 

3.35(15) 

nu 

CsH 2 P0 4 

-260 

- 


0 

- 

HU 

BP0 4 

-0.92(10) 

2.83(12) 

229(9) 

0-56 

5.5(5) 

H] 

BAsCL 

-1.48(15) 

3.64(11) 

181(4) 

0-51 

6.2(7) 

HD 


without any subsequent structural analysis. HO To some extent 
this paucity of crystallographic data is surprising given the in¬ 
tense interest in the high-pressure ferroelectric behaviour of 
this material . 92 In the context of NLC it is perhaps even more 
remarkable because ultrasonic velocity measurements suggest 
a p —» 0 compressibility of —260 TPa -1 along a direction ap¬ 
proximately aligned with the c axis of the monoclinic cell . 87 
The powder diffraction patterns of Ref. 90 do not show any 
obvious evidence of extreme NLC, although re-measurement 
across a larger number of more finely-spaced pressure inter¬ 
vals would help settle the issue definitively. The most likely 
mechanism responsible for NLC in the system (as suggested 
in Ref. [8Tb would involve PO 4 rotation-driven collapse of 
the herringbone hydrogen-bonding network as illustrated in 

Fig-|9fb). 

Because extended structures tend to be less dense than their 
tilted counterparts, more often than not it will be the case 
that pressure-driven activation of tilt systems will favour PLC 
rather than NLC. Cristobalite-like BEO 4 (E = P, As) frame¬ 
works are an interesting case of indirect NLC materials where 
NLC arises because a tilt-driven PLC mechanism results in a 
more rapid compression of the crystal lattice than can sup¬ 
ported by the bulk material stiffness.® We proceed to ex¬ 
plain this mechanism in more detail. Both BPO 4 and BASO 4 
adopt the same tetragonal 74 variant of the cristobalite struc¬ 
ture [Fig. |Tq| . 93 94 The dominant deformation mechanism in¬ 
volves correlated tilting of the BO 4 /EO 4 tetrahedra around the 
tetragonal axis. The tilt angle 0 is directly related to the a lat¬ 
tice parameter via the projection of the mean B-O/E-O bond 
length onto the (a,b) plane r±: S9 



ideal-high- 

(a) cristobalite (b) 


M A 




its 


p = 6.1 GPa 




Fig. 10 (a) The structure of BEO 4 (E = P, As) compounds is related 
to that of cubic j 8 -cristobalite via decoration of the lattice with 
alternating BO 4 and EO 4 tetrahedra and subsequent rotation of all 
tetrahedral units about the tetrad axis by an arbitrary angle 
0 < 0 < 7 t/4 . 93 94 (b) The structures of BPO 4 at high pressures are 
described by increasingly large values of 0. 8 -^ (c) That variations in 
0 dominate the compression mechanism is evident from a 
comparison of the actual variation in the lattice parameter a (open 
squares) with that obtained from the variation in 0 via Eq. (20) 
(solid line). 


ciently large that the inequality 


a 


2 


16 r\ 

tan 2 <j) + 1 ' 


( 20 ) 


Because bending of the B-O-E bonds carries a lower energy 
penalty than compression of the B-0 or E-0 bonds, the be¬ 
haviour of the a lattice parameter on compression is domi¬ 
nated almost exclusively by changes in 0. What this means 
is that K a is effectively a measure only of the bending stiff¬ 
ness of the B-O-E linkages. In contrast, the bulk modulus 
B measures the resistance to compression of the whole oxide 
lattice, which will depend largely on anion repulsion in three 
dimensions. For both BPO 4 and BASO 4 the average value of 
this bulk modulus over the pressure range studied (. B ) is suffi- 


K c = ^-2K a <0 (21) 

holds and NLC is observed along the tetragonal axis [Ta¬ 
ble^ . 89 More remarkable than the magnitude of NLC in these 
systems is the pressure range over which it is observed: Ref. 
^reports K NLC = -0.92(10) and -1.48(15) TPa -1 for E = P 
and As, respectively, over a pressure range of 0-52(5) GPa. 

3.3 Helices 

In many ways the mechanisms covered above—ferroelastic 
instabilities and polyhedral tilting—are similar to those often 
invoked in descriptions of other mechanical anomalies such as 
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Fig. 11 (a) The crystal structure of selenium and tellurium consists 
of a triangular array of trigonal helices (shown in green), (b) The 
structure is completely described by the three parameters r,R, 6 as 
described in the text, (c) Lattice parameter variation under 
hydrostatic pressure for Te, showing NLC parallel to the helix 
axis.®? 103 (d) The dominant compression mechanism involves 
reduction of inter-helix separation R and unwinding of the helices 
(i.e., increasing 0) — it is the latter that gives rise to NLC. 


NTE . 95 Over the next two sections our focus shifts away from 
these dynamical mechanisms to a consideration of topological 
motifs that show a similar predisposition towards NLC. 

Nature herself favours one particular geometric motif— 
namely, the helix—as a mechanism of generating and exploit¬ 
ing NLC in muscular response . 96 Just as NLC materials ex¬ 
pand under hydrostatic pressure, so do they contract under 
negative (i.e., internal) pressure. When filled with fluid, he¬ 
lices exhibit precisely this response: their length decreases as 
the helix cross-section increases. Helical arrays of tendons 
enable muscle-like contraction that is driven by fluid injection 
rather than electrical impulse, and such a mechanism is impli¬ 
cated in the movement of certain types of worms, squid, and 
ancient limbless tetrapods . 96 The NLC community has long 
been aware of the related implication that non-biological NLC 
materials might be exploitable as artificial muscles and actu¬ 
ators if these motifs can be incorporated as part of materials 
design . 97 " 101 

On the atomic scale, there are two remarkably simple chem¬ 
ical systems with helical structures that exhibit NLC under 
hydrostatic pressure. These are the trigonal polymorphs of 
elemental selenium and tellurium . 104 Both structures consist 
of a array of (enantiomorphic) trigonal helices packed on a 
triangular lattice [Fig. ED- Bonding interactions within any 
given helix are much stronger than those between helices so 


an interpretation of the bonding as molecular still has rele¬ 
vance despite the semiconducting properties of both systems. 
The densification required under increasing pressure can then 
be accommodated in two ways: either through compression 
of the weaker bonds between helices—which acts to decrease 
a and leaves c unchanged—or via an increase in the ‘pitch’ 
of each helix, compressing a at the expense of some expan¬ 
sion along c. That both mechanisms operate for Se and Te is 
evident from variable-pressure crystallographic measurements 
[Fig. 0 20 105 j s conventional to reparamaterise the three 
measurable structural parameters— a,c, and u , the single free 
positional variable for Se/Te—in terms of the geometric pa¬ 
rameters 


r 

R 

0 


3 (ua) 2 



[a 2 (l-3u)+ r 2 ] 1 ^ 2 , 
[2cos~ l (3ua/2r)] . 


( 22 ) 

(23) 

(24) 


Here r, R , and 0 correspond, respectively, to the strongly- 
bonded E-E distance, the nearest weakly-bonded E-E dis¬ 
tance, and the E-E-E angle within a strongly bonded helix . 104 
Experimentally one finds that R decreases most rapidly on in¬ 
creasing pressure: this reflects the rod packing compression 
mechanism described first above. In contrast, r is essentially 
constant, meaning that the length of the helical path is essen¬ 
tially incompressible. By themselves, these two observations 
would not imply NLC. Instead the weak NLC behaviour that 
does occur [Table [2) can only be a consequence of the gentle 
increase in bond angle 6 at higher pressures— i.e., the same 
pitch-variation mechanism Nature uses in muscle contraction. 

Remarkably, this NLC mechanism persists even when Se 
chains are incorporated within the cavities of zeolite AIPO 4 - 
5 single crystals E 06 Interpretation of Raman spectroscopy 
and optical absorption measurements, together with ab ini¬ 
tio calculations, converges on the same conclusion that the 
confined chains of Se do actually elongate under hydrostatic 
pressure.EH Indeed, one anticipates that NLC may actually 
be quite a general phenomenon for ordered structures with 
helical motifs. It has been suggested elsewhere that the 
filamentous supramolecular polymer {Au[(C 2 H 5 ) 2 NCS 2 ]} w • 
XCH 2 CI 2 may exhibit negative compressibility via such a 
mechanism ; 23 likewise helical metal-organic frameworks 
such as silver(I) dicyanamide (Refs. 110711081) may be inter- 


Table 3 Compressibilities of helical NLC materials. 


K a 

K c 

Range 

%K 

Refs. 

(TPa -1 ) 

(TPa -1 ) 

(GPa) 

(%) 


Se 12.0(6) 

-2.5(4) 

0-5.2 

1.3(2) 

1201105 

Te 13.6(9) 

-1.8(3) 

0-4 

0.72(10) 

[I 031 03] 
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[NHJ[Zn(HCOO) 3 ] Zn[Au(CN) 2 ] 2 



Fig. 12 Some wine-rack and honeycomb-like topologies known to favour NLC (top) and corresponding chemical systems (bottom). In all 
cases the mechanism responsible for NLC is analogous: densification involves extension of the lattice in the vertical direction. 


esting candidates for further investigation. 


3.4 Molecular Frameworks 


From an engineering perspective—rather than a biological 
one—the two geometric motifs most frequently associated 
with negative compressibility are the wine-rack and honey¬ 
comb networks (which are of course related to each other; 
see Fig. 12). 109 110 Both are characterised by extreme me¬ 
chanical anisotropy and both have the property that their vol¬ 
umes are reduced under uniaxial expansion. Molecular frame¬ 
work chemistry offers an attractive means of designing mate¬ 
rials that contain these same geometric features on the atomic 
scale. 111112 Not only does the chemist have control over net¬ 
work topology, but the use of molecular linkers in the con¬ 
struction of framework materials leads to systems for which 
mechanical deformation mechanisms are dominated by frame¬ 
work flexing. 113 lIt is this feature that enables the mapping be¬ 
tween chemical and engineering systems. 

Within the field of metal-organic framework chemistry, the 
family of materials with a crystal structure most obviously 
resembling the wine-rack topology is probably the so-called 
MIL-53 system. 114 NLC behaviour was predicted for this 
family from both first principles 37 38 and simple empirical 39 
considerations; the first experimental verification of this be¬ 
haviour has only recently appeared in print. 115 The extreme 
magnitude of NLC—Ref. lll51 reports a compressibility value 
of £nlc = — 27 TPa -1 over 0-2 GPa—reflects the shallow en¬ 


ergy potential that accompanies framework deformation for 
these systems and which is implicated in the well-known 
“breathing” effect on guest sorption. 116 Within the field there 
is probably an expectation that such large values are likely a 
general result for molecular frameworks, arising as a conse¬ 
quence of their low-density structures and the low energies of 
the supramolecular interactions affected by framework flex¬ 
ing. 

Our own entry point into this field was via the 
equally remarkable material silver(I) hexacyanocobaltate(III), 
Ag 3 [Co(CN) 6]. 117 Its trigonal structure can be considered a 
three-dimensional wine-rack and is topologically equivalent to 
three interpenetrating a-Po (cubic) nets.L 118 Originally stud¬ 
ied for its bizarre thermal expansion behaviour, the mate¬ 
rial was shown to admit very large strains (^ 10%) even un¬ 
der mild conditions. 19 117 An extremely strong NLC effect 
(^nlc = —76(9) TPa -1 ) was identified in a subsequent neu¬ 
tron scattering study. 19 The mechanism responsible for NLC 
was straightforward enough: densification of the winerack¬ 
like framework proceeded via rapid compression of the a and 
b crystal axes and expansion along c. Despite the magnitude 
of this NLC response, it is unlikely to find widespread appli¬ 
cation since in maximising |£nlc| is framework becomes es¬ 
pecially sensitive to shear instabilities. Indeed a shear-driven 
collapse occurs at p = 0.19 GPa, resulting in rapid densifica¬ 
tion and order-of-magnitude reduction in NLC behaviour. 19 
So while it seems that molecular frameworks can access much 
more extreme compressibilities than the “conventional” solid 
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state materials described in the previous section, there is a 
clear challenge associated with balancing large flexibility with 
a propensity for mechanical instabiliity. 

One strategy for extending the range over which NLC is 
observed is to frustrate collapse by inclusion of counterions 
within the framework cavities. 22 In the case of Ag 3 [Co(CN)6l, 
for example, this can be achieved by substitution of Mn 2+ 
for Co 3+ while balancing charge with extra-framework K + 
ions D 23 The resulting compound KMn[Ag(CN) 2 l 3 remains 
mechanically stable up to at least 2.2 GPa (the highest pres¬ 
sure for which diffraction data have been measured for this 
system), while retaining respectable NLC properties in the 
process: £nlc = —12.0(8) TPa -1 over the range 0-2.2 GPa. 22 
Once again the NLC mechanism involves hinging of the mate¬ 
rials’ wine-rack-like structure. The hypothesis that counterion 
inclusion resulted in “soft-mode frustration” has recently been 
verified directly using a combination of Raman spectroscopy 
and first principles calculations, n 3 and bears similarity to the 
sorption-induced stiffening of metal-organic frameworks doc¬ 
umented elsewhere. 124 

Already some design rules for maximising NLC begin to 
emerge: 

1. Mechanical anisotropy is clearly key: NLC requires PLC 
in an orthogonal direction since volume compressibility 
must remain positive. 

2. Maximising this positive compressibility will likely max¬ 
imise NLC. 

3. Dynamic instabilities will reduce the pressure range over 
which NLC can be observed: using rigid molecular 
linkers and occupying void space with extra-framework 
cations or sorbate molecules may help extend this stabil¬ 
ity range 

4. The network topology should likely be related to the 
wine-rack or honeycomb motifs. 

These design principles were used to identify the first of 
the so-called “ gian t” NLC compounds: zinc dicyanoaurate(I), 
Zn[Au(CN)2]2- The term “giant” demarcates exceptionally 
strong NLC (£nlc < —30 TPa -1 ) that persists over an indus¬ 
trially relevant pressure range (at least 1 GPa). The quartzlike 
structure of this particular material is at once both anisotropic 
and related to the honeycomb net, satisfying respectively the 
first and fourth principles and so explaining the basic driving 
force for NLC. 125 The structure also has very little free vol¬ 
ume because six of these nets interpenetrate one another; 125 
this addresses the third principle and accounts for the extended 
pressure range over which NLC is observed (0-14.2 GPa). But 
what sets Zn[Au(CN)2]2 apart is the strength of NLC is in 
its ambient phase (£nlc = —42(5) TPa -1 over 0-1.8 GPa), 



Fig. 13 Compressibility behaviour of a variety of molecular 
frameworks. The width of each bar is related to the pressure range 
over which compressibility is measured, and the height of each bar 
corresponds to the maximum (red) and minimum (blue) average 
compressibility values over that range. Consequently the area of 
each bar is a measure of %K- those materials with the greatest 
propensity to expand under pressure are those for which the area of 
the blue bar is largest. 


which arises from the extreme compressibility of supramolec- 
ular helical “springs” allowing especially strong PLC in one 
set of directions. Flexing of the honeycomb-like framework 
translates this exceptionally strong PLC into an equally re¬ 
markable NLC effect along the c crystal axis. It is this use 
of a supramolecular helix that addresses the second design 
principle: just as a spring is more compressible than the steel 
from which it is made, so too are the linear compressibilities 
of Zn[Au(CN)2l2 more extreme than would otherwise be ex¬ 
pected. 


NLC is increasingly frequently identified in a wide range 
of MOFs and molecular crystals beyond the various systems 
discussed above. Two relevant MOF examples are ammo¬ 
nium zinc(II) formate, [NH 4 ][Zn(HCOO) 3 ]P and the zinc 
alkyl gate (ZAG) family. 120 The mechanism responsible for 
NLC in the former is essentially the same as that discussed for 
KMn[Ag(CN) 2 ] 3 ; all that differs between the two is their net¬ 
work topology (cag vs a-Po). 7 lln contrast, the “wine-rack” 
NLC mechanism originally proposed for ZAGs (Ref. 1201 
was later shown to be incorrect on the basis of quantum me¬ 
chanical calculations. EH Instead it seems that NLC in this 
system is discontinuous and is driven instead by pressure- 
driven proton redistribution—a new mechanism altogether. 
The NLC behaviour of these MOFs is compared with that 
of other molecular frameworks in Table [4] and Fig. 13 In 
terms of molecular crystals, strong NLC effects can occur 
in situations where packing arrangements mimic the topo- 
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Table 4 Compressibilities of NLC molecular frameworks and molecular solids. 



Ki (TPa -1 ) 

K 2 (TPa -1 ) 

K 3 (TPa" 1 ) 

Bo (GPa) 

Range (GPa) 

Z* (%) 

Ref(s). 

Ag(mim) 

-4.32(10) 

25.8(10) 

55(4) 

6.0(16) 

0-1.0 

0.432(10) 

EH 

KMn[Ag(CN) 2 ] 3 

-12.0(8) 

33.2(13) 

33.2(13) 

12.7(11) 

0-2.2 

2.64(18) 

22119 

Ag 3 [Co(CN) 6 ]-I 

-76(9) 

115(8) 

115(8) 

6.5(3) 

0-0.19 

1.44(17) 

EH 

Ag 3 [Co(CN) 6 ]-II 

-5.3(3) 

9.6(5) 

15.2(9) 

11.8(7) 

0.19-7.65 

4.0(2) 

EH 

Zn[AufCN) 2 ] 2 -I 

-42(5) 

52(6) 

52(6) 

16.7(16) 

0-1.8 

7.6(9) 

EH 

Zn[Au(CN) 2 ] 2 -II 

-6(3) 

16(5) 

16(5) 

27(3) 

1.8-14.2 

7(4) 

EH 

[NH 4 ][Zn(HCOO) 3 ] 

-1.8(8) 

15.8(9) 

15.8(9) 

32.8(16) 

0-0.93 

0.17(7) 

m 

ZAG-4 

-2.6(15) 

7.9(5) 

29(3) 

11.66 

1.65-5.69 

1.1(6) 

E 20 I 

[Fe(dpp) 2 (NCS) 2 ]-py 

-10(2) 

12(3) 

53(4) 

12.9(6) 

0-2.48 

2.5(5) 

11211 

CH 3 0HH 2 0 

-2.7(18) 

31.9(4) 

108.0(9) 

3.79(6) 

0-0.6 

0.16(11) 

E 22 I 

[(C6F5Au)2(jU-l,4-diisocyanobenzene)] 

-13(3) 

29(3) 

31.3(4) 

7.5(7) 

0-2.42 

3.1(7) 



logical motifs known to favour NLC in framework struc¬ 
tures. 127 So, for example, the molecular packing arrangements 
in systems as chemically diverse as methanol monohydrate, 
[Fe(dpp) 2 (NCS) 2 ]-py (dpp = dipyridoP^-^^'-cjphenazine, 
py = pyridine) and [(C6F5Au)2(/i-l,4-diisocyanobenzene)] 
are all related to the same wine-rack topology. Mod¬ 
erately strong NLC is observed in each case: £nlc = 
—2.6(3), —10.3(20), and —4.16 TPa -1 , for each of these three 
examples in turn. 9 121122 

4 Discussion and Future Directions 

At the time of Baughman’s original review (Ref. |2]), NLC ap¬ 
peared to be a somewhat esoteric phenomenon that occurred in 
only a handful of peculiar systems. If the enumeration of vari¬ 
ous NLC materials and mechanisms given in our review leads 
to any one particular conclusion it is surely that NLC is rather 
more commonplace than might originally have been expected. 
From an experimental viewpoint, characterisation of NLC re¬ 
mains somewhat of a niche capability; however, the number 
of research groups with expertise in high-pressure crystallo¬ 
graphic measurements is rapidly growing and we anticipate 
that NLC will be increasingly frequently observed as the scope 
of materials studied at high pressures diversifies. This is likely 
to be particularly true for molecular crystals, where the empir¬ 
ical predisposition towards herringbone-type packing arrange¬ 
ments will intrinsically favour NLC because of its relationship 
to the winerack topology. 

But perhaps the most profound recent advance likely to 
influence our understanding of the true breadth of NLC be¬ 
haviour is the relative ease and reliability with which first 
principles calculations can now determine elastic properties. 32 
This allows (relatively) rapid screening of entire classes of 
materials and pre-selection of the most interesting candidates 
for subsequent experimental investigation. We expect that the 
task of calculating elastic tensors for all known MOF struc¬ 
tures, for example—once considered an inconceivably diffi¬ 


cult process—will not remain computationally intractable for 
very much longer. 126 

If NLC is likely to be discovered in a large number of 
new systems over the coming years, then what are the cri¬ 
teria for interesting and useful behaviour? Magnitude and 
range of NLC are two obvious metrics—hence the otherwise- 
questionable value of labels such as “giant” NLC. 40 But one 
other aspect of NLC behaviour deserves brief discussion: 
namely the relationship between NLC and crystal symmetry. 
Low symmetry intrinsically favours anisotropy—and hence 
NLC—since there are fewer constraints on the elastic stiff¬ 
ness tensor C. While NLC is likely to be more prevalent in 
low-symmetry materials, the direction along which NLC oc¬ 
curs is less likely to have a fixed relationship to the crystal axes 
(hence bulk sample morphology) and moreover this direction 
itself will vary with pressure. 15 This matters because practical 
implementation of NLC requires careful material alignment. 
So, in general, we anticipate that NLC behaviour in uniaxial 
or orthorhombic systems will find greatest application because 
in these cases macroscopic alignment of the NLC axes can be 
assured. 

Having documented NLC in large variety of different chem¬ 
ical families, we sought to establish whether any universal 
trends in behaviour might emerge. Our approach is to compare 
for these different materials the relationship between magni¬ 
tude of NLC effect and pressure range over which the phe¬ 
nomenon is observed. Our analysis is summarised in the 
Ashby-type diagram given in Fig. [14] The diagonals in this 
representation correspond to points of constant %k- What 
becomes obvious is that most NLC materials are distributed 
around the diagonal corresponding to %k = 1 %—this value is 
probably special only in the sense that it is large enough for 
NLC to have been noticed experimentally, but is not so large 
as to be extremely rare. Along this diagonal, those materials 
for which the mechanically-relevant bonding interactions are 
strong—the rutiles and framework silicates, for example—are 
clustered in a region where NLC pressure range is large but 
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Fig. 14 Ashby-type diagram relating magnitude of NLC behaviour 
to pressure range over which it is observed for those materials 
included in this review. Diagonals correspond to points of constant 
compressibility capacity. 


^nlc itself is not so very extreme. In contrast, the flexible, 
open structures of MOFs and molecular frameworks cluster in 
a region where NLC is very strong but persists only over small 
pressure ranges. In this way, eventual materials selection will 
be informed by whether range or magnitude is more critical 
for a given particular application of NLC. In terms of new ma¬ 
terials discovery, the most attractive region of the Ashby plot 
is of course the (currently empty) top-right corner, where both 
metrics are maximised. Somehow the materials that exist in 
this region—if indeed they exist at all—must balance the weak 
interactions needed to produce strong NLC with the structural 
integrity required to avoid collapse at high pressures. 

Our focus in this review has been almost entirely on ma¬ 
terials which exhibit negative compressibility along just one 
principal axis. Negative area compressibility (NAC) is a nat¬ 
ural extension and indeed modest NAC effects have been ob¬ 
served or predicted to occur in a handful of layered materi¬ 
als: silver(I) tricyanomethanide,® 8 sodium vanadate, 129 133 
and TlGaSe 2 (Ref. 134) (Table[5j we note that reports of NAC 
in PbTi 03 remain contentious and so are not included 135 137 ). 
In each of these cases, the mechanism responsible forJS1AC 
is related to the so-called Lifshitz mode [Fig. [l5j. 138 An 
alternate mechanism—in which specific geometries of the 
various wine-rack and honeycomb-like topologies identified 
above couple densification to area expansion—has recently 
been proposed;^ the extension to specific helical geometries 
also follows. Whatever the mechanism, it is reasonable to ex- 


Table 5 Compressibilities of NAC materials 



An.AC 

(TPa -1 ) 

KpLC 

(TPa -1 ) 

Range 

(GPa) 

Ref(s). 

Ag[C(CN) 3 ] 

-7.5(8) 

66(20) 

0-0.615(6) 

11281 

NaV 2 0 5 

-3.128 

1.67 

2-10.0 

nmm 

TlGaSe2 

-4.99 

33.28 

- 

[1341 


pect that NAC is unlikely ever to be as strong an effect as 
NLC. This is because the mechanical stability criterion of pos¬ 
itive volume compressibility implies that the PLC effect along 
the axis perpendicular to the plane of negative compressibility 
must be at least twice as large as £nac- In other words, NAC 
can only ever be half as strong as PLC. 

While continuous volume compressibility must be positive, 
there is no formal thermodynamic requirement forbidding vol¬ 
ume increase across a pressure-induced phase transition. To 
the best of our knowledge, the only ‘realisation’ of this bizarre 
phenomenon is a theoretical study of a fictitious multicompo¬ 
nent metamaterial assembled from cleverly-chosen pairwise 
potentials. 140 Translating this study into real materials that 
exhibit a negative compressibility transition is difficult (per¬ 
haps impossible) for three main reasons. First, the expanded 
phase is only stable at high pressures in the 0 K limit and is 
metastable at finite temperature. Second, the barrier to relax¬ 
ation of this metastable state scales inversely with system size. 
For atomic-scale realisations— i.e. materials as discussed in 
this review—the barrier is unlikely to be sufficiently high to 
prevent rapid relaxation to the thermodynamic (compressed) 
state at high pressures. 140 And, third, the pairwise potentials 
needed to produce the effect are as unusual as they are spe¬ 
cific. Nevertheless, were it possible to realise even a discontin¬ 
uous negative compressibility transition in a chemical system 
then “continuous” negative volume compressibility behaviour 
might be feasible through inhomogeneous chemical doping to 
“smear out” the transition pressure—a trick used elsewhere to 
convert discontinuous thermal volume collapse to “colossal” 



P 



Fig. 15 Lifshitz mechanism for NAC. 138 Densification of layered 
materials usually proceeds via collapse in the stacking direction, 
which in turn results in an expansion in the two perpendicular 
directions ( i.e ., within the layer). 
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NTE.ED 

Whatever the particular manifestation of negative com¬ 
pressibility, what is clear is that its very existence is symp¬ 
tomatic of anomalies in the whole elastic tensor. One obvious 
area for development in the field involves the use of elastic 
property measurement techniques—historically applied in the 
fields of materials science, engineering, and metallurgy—as 
a means of exploring related elastic anomalies in the various 
chemical systems covered in this review. One expects that 
the general trends identified here will resurface irrespective 
of the particular measurement: strongly-bound materials will 
be more resistant to applied stress, and open frameworks will 
give rise to larger magnitude elastic responses. The develop¬ 
ment of new classes of negative Poisson ratio materials, for 
example, would be one such avenue of further research. 142 

From a materials chemistry perspective, however, perhaps 
the most exciting direction to be explored will be the coupling 
of anomalous elastic behaviour with other materials proper¬ 
ties. The potential application of NLC in pressure sensing de¬ 
vices, for example, relies on the effect of simultaneous elonga¬ 
tion and densification on optical behaviour. 2 Orbital overlap in 
superconductors and magnetic materials—itself strongly sen¬ 
sitive to the distances between atoms—may also be tuneable 
in unexpected ways by exploiting the counterintuitive pressure 
dependencies attainable of NLC materials. 143 Thinking be¬ 
yond NLC itself, the coupling of extreme elastic anisotropy (of 
which NLC is of course a simple consequence) with macro¬ 
scopic dipole formation may lead to exciting new classes 
of piezoelectrics, ferroelectrics, and pyroelectrics. Whatever 
the future holds, the lessons learned from the study of NLC 
materials—namely the different ways in which chemical, geo¬ 
metric, and topological motifs influence elastic anomlies—are 
likely to play a crucial role in functional materials design. 
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